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SECTION  1 


INTRODUCTION 


The  purpose  of  this  document  is  to  provide  a standard  algebraic 
model  for  a continuous-wave  <CW)  high-energy-laser  (HEL)  beam  on  target. 
The  been  is  degraded  by  thermal  blooming  as  well  as  numerous  other 
effects.  The  algebraic  model  is  intended  for  systems  analysis  exer- 
cises, where  the  numerous  parameters  are  to  be  explored,  making  more 
detailed  computer  simulation  of  atmospheric  propagation  impractical. 

It  is  an  outgrowth  of  efforts  to  fit  scaling  laws  to  a large  body  of 
data  generated  by  detailed  atmospheric  propagation  simulations  using 
a finite-difference  wave-optics  code.* 

Any  such  wave-optics  code  provides  a detailed  intensity  profile 
on  target,  from  which  several  different  summary  characteristics  can  be 
extracted.  These  include,  for  instance,  peak  intensity,  line-of-sight 
beam  dispersion,  and  beam  area.  The  last  of  these  further  admits 
several  definitions,  including  area  to  some  percentage  of  total  power, 
and  area  defined  by  a ratio  of  squared  integral  of  intensity  to  integral 
of  intensity  squared  (suggested  independently  by  Lincoln  Laboratory 
and  Draper  Laboratory  researchers).  If  beams  on  target  were  Gaussian 
in  shape,  all  such  characteristics  would  convey  equivalent  information. 
The  simplified  algebraic  model  assumes  that  this  is  nearly  the  case, 
and  speaks  nominally  of  peak  intensity  on  target. 

Regression  analysis  of  the  results  has  shown  that  peak  intensity 
can  be  correlated  with  an  integral  which  represents  the  accumulation 
along  range  (starting  from  the  center  of  the  aperture)  of  phase  pertur- 
bation due  to  heating  in  a beam  with  absorption,  scattering,  convectiva 
clearing,  and  focusing.  Figures  1,  2,  and  3 show  the  tightness  of  the 
correlation  obtained  for  three  different  beam  shapes.  The  ordinate  is 
the  phase  integral  ¥h  and  the  abscissa  is  the  ratio  R * (Iy  - I) /I, 
where  I is  peak  intensity  and  ly  is  unbloomed  peak  intensity. 


Data  provided  by  D.  Cordray  of  Naval  Research  Lab  (NRL) . 


1 


Any  one  of  many  numerically  aimtlar  functional  forms  could  be  fit 
to  the  curves  defined  by  the  correlation  data.  Thir  document  selects 
one  that  has  been  found  most  useful  because  of  its  algebraic  simplicity. 
Section  2 shows  how  the  simple  algebraic  form  leads  to  a number  of 
universal  relationship*  that  are  independent  of  almost  all  physical 
details  of  the  propagation  process.  These  relationships  describe  the 
variation  of  peak  intensity  as  a function  of  pover  alone,  with  all 
other  variables  held  constant,  important  relationships  are  shown  to 
depend  only  on  beam  shape.  This  conclusion  has  not  previously  been 
evident  from  other  scaling  laws. 

Section  3 discusses  the  parameters  appearing  in  Section  2,  showing 
how  each  depends  on  actual  physical  variables  that  describe  the  engage- 
ment, the  laser,  and  the  atmosphere.  Section  4 provides  a technique 
for  accurately  evaluating  the  phase  integral  along  range  required  in 
Section  3.  For  several  common  beam  shapes.  Section  5 discusses  numerical 
values  of  regression  parameters  defined  in  Section  3.  Section  6 pro- 
vides a concise  summary  of  all  formulas  required  for  systems  analysis 
exercises. 


Figure  1.  Correlation  for  infinite 
Gaussian  beam. 
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Figure  3.  Correlation  for  uniform  beam. 


SECTION  2 


INTENSITY  DEGRADATION  AS  A 
FUNCTION  OF  POWER 


For  a class  of  problems  differing  only  in  time-average  power  (P) , 
the  ratio  of  bloomed  to  unbloomed  intensity  focused  on  target  can  be 
modeled  as 


I (Pi 

V? T 


2 

In  this  expression,  oT  is  the  1/e  line-of -sight  beam  dispersion  in 

radians  due  to  linear  effects,  including  diffraction,  beam  quality, 

turbulence,  and  jitter.  The  latter  two  effects  separate  irfco  "high" 

and  "low"  frequencies,  vhicn,  respectively,  do  and  do  not  impact 

2 

blooming.  The  latter  is  represented  in  the  term  oQ.  To  some  extant, 

the  form  of  X(P)/Xg(P)  expresses  the  familiar  idea  of  "rss-ing"  (root- 

sum-squaring)  . For  a Gaussian  beam,  dispersive  effects  combine  by 

2 2 

summing  variances,  and  the  denominator  (o^  + afi)  resembles  such  a sum. 

2 

The  blooming  term  (oQ)  depend j on  power  (P)  in  a way  that  can  be 
modeled  by  a variety  of  functional  forms.  The  choice  of  form  is  a 
tradeoff  between  simplicity  and  range  of  validity.  In  this  report, 
we  use  the  simplest  form  Known  to  be  valid  for  power  levels  of  practical 
interest.  This  form  is  nominally 


o 


2 

B 


C P 


The  use  of  Pa  with  a > 1 allows  reproduction  of  a well-known 
physical  phenomenons  there  exists  a critical  power  Pc  such  that  the 


* 


This  type  of  formula  has  also  been  suggested  and  used  by  F.  Gut hard t 
and  J.  Wallace.  In  particular,  Wallace  suggested  a - 3/2. 


1 

intensity  on  target  is  maximum.  That  Pc  can  be  related  to  o^,  CB, 
and  a as  follows.  The  unbloomed  intensity  on  target  is  proportional 
to  P,  so  that  actual  intensity  is  proportional  to  P[1 (P)/Zn(p) ] , or 
Po£/(o‘  + CgF  ).  Differentiating  with  respect  to  P and  setting  the 
derivative  to  sero  gives 


ll/a 


(a  - ire. 


From  the  value  of  Pc,  the  ratio  X(Pc)/zu(Pc)  is  readily  found 

to  be 


l(Pc) 

W 


a - 1 


a 


This  result  is  interesting  because  it  depends  only  on  the  parameter  a, 
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expressive  r.f  beam  shape  alone,  and  not  on  or  CB,  which  contain  much 
of  the  physics  of  the  problem.  Zt  therefore  presents  a physical 
phenomenon  that  is  essentially  separable  from  other  physical  phenomena 
in  the  overall  propagation  process. 

The  result  concerning  I(PC)/Z0<PC)  indicates  what  to  axpect  from 
adaptive  phase  correction  for  thermal  blooming.  Since  phase  correction 
does  not  change  beam  shape,  it  will  move  the  whole  curve  of  intensity* 
out  versus  power- in  in  such  a way  that  the  new  peak,  the  old  peak,  and 
the  origin  lie  on  a straight  line,  as  illustrated  in  Figure  4. 


Clearly,  Pc  marks  the  upper  limit  of  power  levels  having  practical 
interest.  In  fact,  operation  well  below  ?c  may  be  of  interest,  so  let 
us  consider  P at  aoeie  fraction  of  ?ci  P • Pc/b.  Then 


which  simplifies  to 


* W* 


ab 


a-1 


b (a  - 1)  + 


1 


2 

This  result,  too,  is  independent  of  and  CB,  and  therefore  is 
independent  of  the  phenomena  controlling  them. 
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Figure  4.  Adaptive  correction  of  phase  without 
change  of  beam  shape. 

The  formula  describing  operation  below  Pc  is  plotted  in  Figure  5 
for  several  typical  values  of  a.  For  all  cases,  there  is  little  point 
in  operating  at  power  levels  above  Pc/2,  because  more  than  85»  of  the 
limiting  peak  intensity  is  already  available  at  Vq/2. 

The  formulas  for  I(P<;)/1U(PC)  and  1 (Fc/b)/X  (Pc)  can  be  combined 
to  relate  I (P^yb)  to  Iy  (Pc) . The  result  is 


£ 


(a  - l)ba‘X 
b*(a  - 1)  + 1 


This  result,  too,  is  independent  of  the 
most  of  the  physics  of  the  prohlem.  Zt 
is  specified,  I for  any  other  condition 


2 

parameters  oL  and  CB  that  carry 
means  essentially  that  if  IU(PC) 
is  implicitly  specified. 


j 
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The  nodal  presented  in  this  section  is  most  accurate  in  the  regime 

of  physical  interest  for  systems  modeling,  namely,  below  the  sero-wander 

value  of  Pc<  This  is  the  case  because  below  this  power  there  is  relatively 

little  blooming,  with  little  attendant  beam  distortion  and  no  beam 

breakup.  There  is  then  a one-to-one  Gaussian- like  correspondence  between 

2 

intensity  I and  the  various  disf"  rsivo  o terms.  Operation  far  above 
such  a power  level  requires  a more  complex  model,  involving  linear  and 

j (1)* 

quadratic  laser  power  terms  in  o^,  presented  by  Breaux. 

The  ratio-type  results  of  this  section  are  of  a very  general 

nature,  and  are  valid  regardless  of  the  methodology  used  to  evaluate  a, 

2 

; oL,  CB,  or  Ig(Pc).  Specific  techniques  for  evaluating  these  are 

provided  in  the  Sections  3 through  5. 

S t ^ 

F Superscript  numerals  refer  to  similarly  numbered  items  in  the  List 

[ of  References. 
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SECTION  3 


PHYSICAL  MODEL 


The  model  for  degradation  of  intensity  focused  on  target  as  a 

2 

function  of  power  has  only  three  parameters  (a,  o^,  and  CB) , and  requires 
only  I0(PC)  to  specify  intensity  under  any  conditions.  A small  amount 
of  data  for  a given  class  of  problems  differing  only  in  P will  suffice 
to  determine  all  the  required  quantities.  But  extrapolating  the  re- 
sults to  any  other  class  of  problems  requires  some  kind  of  physical 
model.  The  exponent  a can  be  assumed  to  be  independent  of  many  physical 
variables,  so  the  more  pressing  problems  are  to  admit  defocus,  and  to 
model  c£,  CB,  and  Iy (Pc) . 

The  physical  variables  affecting  HEL  propagation  are  designated 
in  Table  1.  In  cases  where  the  propagation  path  extends  through  signif- 
icantly different  altitudes,  path  averaging  of  parameters  is  required. 
(See  Reference  1,  p.  59  for  applicable  techniques  regarding  c\) 

These  parameters  can  be  sunnarized  in  terms  of  four  dimensionless 
numbers: 


Absorption  number 

* 

Slew  number  Ng 


Fresnel  number  Np 

*+ 

Distortion  number  ND 


* 

Defined  appropriate  to  this  model. 

* Subscript  » me'  ns  unperturbed  natural  value. 
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Table  1.  Physical  variables  affecting 
HEL  propagation. 


Parameter 

Variable 

Units 

Engagement 

Target  range 

Zt 

meters 

Transverse  target 
velocity 

Ut 

meters/second 

Laser 

Aperture  radius 

meters 

Time  average  power 

p 

Joules/second 

Wavelength 

X 

meters 

Beam  quality 

M 

no  dimension 

Beam  shape 

no  dimension,  not  a scalar 

Atmosphere 

Wind  velocity 

°w 

meters/ second 

Turbulence 

meters”2^3 

Absorption 

a 

meters  3 

Scattering 

a 

meters”3 

Index  of  refraction 

n 

no  dimension 

Temperature 

T 

degrees 

Refraction  gradient 

3n/3T 

degrees”3 

Density 

P 

kilogram-meters-3 

Heat  capacity 

CP 

Joules/degree-kilogram 

There  are  also  effective  bean  qualities  Nq  and  N£,  which  repre- 
sent the  actual  sheading  of  the  beam  due  to  various  effects.  These 
effects  include  beam  shape  and  the  blooming  itself,  and  so  and 
Nq  must  be  defined  later.  A procedure  is  provided  in  Section  6. 

2 

The  a L is  1/e  radial  beam  dispersion  due  to  linear  effects.  It 
includes  a number  of  contributors,  distinguished  by  different  subscripts: 
diffraction  and  beam  quality  (D) . high-frequency  beam  distortion  and 
motion  due  to  turbulence  (T)  and  jitter  (J) , low-frequency  beam  wander 
(W)  due  to  turbulence,  jitter,  and  pointer-tracker  effects.  Thus,  we 
have 


o 


2 

L 


oj?  + c*  + 

D T 


The  formula  for  the  diffraction  contributor  is 


o 


2 

D 


2 


where  m*  is  characteristic  of  beam  shape.  Requiring  aQ  to  be  the 

63%  beam  radius  on  target  makes  m'  equal  to  2/ir  for  an  infinite 

2 

Gaussian  beam,  equal  to  0.9166  for  a 1/e  truncated  Gaussian  beam,  and 
equal  to  0.9202  for  a uniform  beam  (see  Reference  2). 

The  turbulence  contributor  refers  specifically  to  the  short-term 
high- frequency  part  of  total  turbulence.  It  has  been  found  by  Breaux  ^ 
that  for  a large  variety  of  beam  shapes,  total  turbulence  is  well 
represented  by 


where  rQ  is  the  Fried  coherence  length  for  wave  number  k ■ 2ir/X 
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and  De  is  an  effective  aperture  size,  appropriate  to  the  beam  shape.  For 
infinite  Gaussian,  truncated  Gaussian,  or  uniform  beams,  respectively, 
it  is 


aRm'  3'74'  4R» 


The  short-term  part  (oT)  is  smaller,  varying  between  the  values 


0.182 


for  De/rQ  < 3 and 


2 

/OfA 

Y2  / 

D 

hr) 

G 

')  - U l8( 

V 

for  De/rQ  > 3. 

To  ensure  model  validity,  the  high-frequency  turbulence  should  be 

limited  to  values  small  enough  to  cause  no  speckling,  say  a < 2 on.  The 

2 2 1 2U 
difference  (oTT  - oT)  may  or  may  not  appear  in  the  wander  term  a^, 

depending  on  the  particular  hardware  implementation  being  modeled. 

Additional  jitter  and  wander  contributors  arise  from  the  particular 

hardware  considered,  and  must  be  set  by  the  user. 

The  dispersion  due  to  blooming  is  to  be  combined  with  the  linear 

2 2 2 2 

dispersions  simplified  by  cD,  aT,  aJf  and  a w.  We  proceed  now  to 
consider  the  blooming  term 


a 


2 

B 


CBP 


a 


2 

The  physical  phenomenon  that  causes  is  accumulation  along  range  of  phase 

due  to  heating,  which  increases  with  power  (P).  and  decreases,  to  some 

2 

extent,  with  that  portion  of  that  actually  experiences  the  blooming. 
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2 2 

namely  o£  - aw.  It  has  been  found  possible  to  define  a measure  Vh  for 
heating  phase  such  that  the  following  expanded  expression  trail  represents 


°l  - wl-0 i><va 


where  Cg  is  a dimensionless  coefficient  that  depends  on  beam  shape. 
Clearly 

s - ci<°l  - °£>(r) 

2 

The  form  of  that  makes  the  above  representation  of  possible 
is  proportional  to  NpNp/N^.  The  phase  integral  (¥h)  also  has  an 
additional  factor  to  make  it  into  an  integral  along  range  that  takes 
account  of  extinction,  clearing,  and  focusing  of  the  beam.  Its  evalua- 
tion is  the  subject  of  Section  4. 

Next,  I0(PC)  is  easily  estimated  by  considering  total  extinction 
(e  *=  a + a),  beam  shape  and  quality,  and  the  spreading  due  to  linear 
effects.  The  estimate  is 


W 


Pc  ®xp(-e&t> 
2 2 

2*Z£o£B 


Here,  B is  a beam-shape  factor,  which  can  be  roughly  interpreted  as  the 
ratio  of  average  intensity  to  "peak"  intensity  in  the  focal  plane.  For 
a Gaussian  beam,  it  is  exactly  1/2,  and  for  any  realistic  beam  shape, 
it  is  quite  close  to  1/2. 


fc 


SECTION  4 

PHASE  INTEGRAL 


In  Section  3,  the  problem  of  evaluating  was  reduced  to  the 

problem  of  evaluating  the  phase  integral  T.  . Many  integration  pro- 

(145)  n 

cedures  have  been  investigated,  ' and  most  can  be  made  numerically 
adequate.  The  following  is  straightforward  to  explain.  Let 

1 

’h  ■ / V*>  dI 

o 


where 


z 


V‘> 


NpNp  exp(-eZ)  J 

1 1 

nq 

Lv  J 

3“ 

M 

r2V 

1 

i 

i 


I 


t 


i 


Here,  exp(-EZ)  represents  extinction.  The  variable  Rg(z)  represents 
spot  radius  at  s,  approximated  by  rss-ing  focus  and  diffraction  effects. 


Ra(s) 


1/7 


The  term  29^/V^  represents  the  clearing  time  at  the  aperture,  and  pro- 
vides a normalisation  for  the  time  integration.  P/*R^  provides  a 
normalisation  for  the  intensity,  which  depends  on  time  (t)  through  xQ 
and  Yg,  which  are  both  0 for  t - 0.  For  a beam  slewed  in  the  x 
direction 


j 

\ 


yQ(s,t)  - 0 
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but  XQ<i,t)  varies  with  the  local  clearing  velocity.  Assuming  xQ,  yQ 
nondimensionalised  by  aperture  radius 

-tUJs) 

x0(z'tJ  " 1.1 7T 

with  Vm(x)  representing  that  velocity,  approximated  by 

U„(*)  - UW<1  + sNg) 

As  an  example,  consider  a Gaussian  beam  with  amplitude 


A0 ^ 0 ^ 


exp(-xj  - Vq) 


For  the  Gaussian  profile,  the  normalized  intensity 

1 [x0 (2,t) » y0(2'fc)] 


■ 2 exp 


The  time  integral  of  normalised  intensity  is 


L2vw/ 


2H,(Z)1/Uw^  [i  (l\ 

\2/ 


Absorbing  the  /ir/2  (1/2)  in  the  overall  relating  to  beam  shape,  and 
substituting  for  U^fs)  and  Rs(z),  leaves  the  required  z integral  to 
be  just 


V 


h 


/ 


NDNF 


N£(!  + ZNS) 


exp(-eZ) 


The  integral  does  not  appear  to  be  amenable  to  direct  analytic 
evaluation,  so  the  options  are  numerical  integration  or  analytic 
approximation.  Typically,  numerical  integration  is  difficult  because 


14 


aaall  steps  are  required  to  handle  the  rapid  variation  of  the  integrand. 
Therefore,  analytic  approximation  ia  preferred.  Probably  the  feweat 
approximations  are  required  if  integration  by  parta  ia  employed,  in 
that  caae  we  find 

1 

*h  - Jaxp(-eZ) Y^]1  + exp(-eZ)Y£  da 

0 0 

The  indefinite  integral 


ia  available  in  an  integral  table  by  Klerer  and  Groaeman. ^ Naturally, 
the  expreaaion  containa  an  integration  conatant.  In  principle,  thia 
conatant  needa  to  be  reaet  in  auch  a way  as  to  nullify  the  remainder 
term 

1 

clt  / exp(-eZ) dz 
0 

Thia  can  be  accompliahed  at  leaat  approximately  by  using  Y^(i)  - V^(l/2) 
in  place  of  Y^(z).  Thia  makea  the  expreaaion  for  Yh  quite  complicat'*^. 
Utility  ia  greatly  served  by  reducing  it  to  the  caae  of  small  attenu 
tion  ezt,  where  the  integration  constant  has  no  impact,  and  Y^  reduces 
to 


Tha  condition  for  thia  ra«ult  to  ba  valid  can  ba  atatad  as 


SECTION  5 


CORRELATION  DATA 


Extensive  simulations  of  atmospheric  propagation  have  been  per- 
formed at  NRL  using  a finite-difference  wave-optics  code,  and  these 
provide  the  data  base  that  establishes  a strong  correlation  of  the 
form  suggested  in  this  report  relating  intensity  on  target  to  the 
parameter  There  is  very  little  scatter,  and  that  which  does  exist 

may  be  attributed  to  one  of  two  factors: 

(1)  Earlier  correlation  studies  by  Seiders^  established  that 
the  spot  site  in  the  integrand  denominator  of  should 
not  be  just  the  ideal  diffraction-limited  value,  mainly 
because  of  the  iterative  effect  of  blooming  upon  itself. 

To  account  somewhat  tor  this  in  a way  that  maintains 
reasonable  simplicity,  a free  constant  (m)  scaling  the 
diffraction  spot  was  introduced  and  evaluated  by  minimising 
residuals.  Actually,  such  scaling  must  depend  to  some 
extent  or  physical  variables,  especially  power.  The  idea 
ci  iteration  to  adnit  variable  m has  been  tried,  but  so 
far  has  not  been  sufficiently  successful  to  justify  the 
effort. 

(2)  The  raw  data  generated  were  not  all  of  the  same  form.  Data 

for  the  infinite  Gauasian  beam  comprisad  paaK  intensities, 

2 

whereas  those  for  the  truncated  Gaussian  comprisad  1/a  area, 
and  '-.hose  for  the  uniform  beam  comprised  tha  more  complex 
area  measure*  of  the  form 


[/ [ d>  dyj2 


f 


X dx  dy 


The  latter  functional  (A)  is  currently  thought  to  lead  to 
the  tightest  correlations. 
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2 

The  result-  of  analysing  three  beam  shapes — infinite  Gaussian,  1/e 
truncated  Gaussian,  and  uniform  beam— -are  presented  in  Tables  2 through 
4 and  plotted  in  Figures  6 through  8.  Supposing  that  NQ  and  Np  are 
fixed  and  o*  equals  sero,  let  us  compare  the  three  beam  types.  We 
recall  that  Cfi  is  proportional  to  and  that 


P 


C 


Ti~-"TO 


l/a 


and 


iirc'  a - 1 

The  results  tend  to  follow  intuitively  understandable  patterns.  First, 
the  parameter  m affects  the  tightness  of  the  correlation  much  less  in 
tha  case  of  the  Gaussian  beams  than  in  the  case  of  the  uniform  beam. 

This  is  to  be  expected  because  blooming  in  the  Gaussian  case  is  driven 
mainly  by  gradients  near  the  aperture,  whereas  for  the  uniform  boam  such 
gradients  are  minimised  and  blooming  is  driven  by  beam  shape  nearer 
the  focus,  which  is  described  by  a.  Secondly,  the  uniform  be/m  has 
the  smallest  a and  C^/o^,  resulting  in  by  far  the  largest  P^.  The 
large  Pc  is  to  be  expected  because,  with  less  gradient  near  the  aperture, 
there  is  less  lensing  effect  .even  at  high  power.  The  small  a means, 
however,  that  the  Pc  region  is  very  broad  and  I (Pc)/Iy (fc)  is  very  entail. 
The  two  Gaussian  beams  have  anomalously  different  C£,  a,  and  m,  but 
they  turn  out  to  have  similar  small  P^  values.  But  because  a is  larger 
for  these  cases,  I(Pc)/Ig(Pc)  is  somovhat  larger. 

The  actual  peak  intensity  on  target  at  Pc  depends  on  Pc»  on 
I(PC)/IU(PC)  - (a-l)/a,  and  on  the  diffraction-limited  spot  radius  r( 
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The  spot  radius  for  the  three  bean  shapes,  respectively,  is  proportional 
to  2/ir,  0.3166,  and  0.9202.  with  the  result  that  I(PC)  for  the  truncated 
Gaussian  is  only  half  that  of  the  infinite  Gaussian.  The  value  for  the 
uniform  beam  is  nearly  three  times  what  it  is  for  the  truncated  Gaussian, 
but  at  the  price  of  nearly  six  times  the  power  (which  may  nevertheless  be 
acceptable) . 


Table  3.  1/e  truncated  Gaussian t I(P)/Iy(P) 


(1  + 


rma 

Error 


Maximum 

Error 


Choaen  j 


0.170 

0.169 


2 

8 


0.981 

0.924 


0.018600 

0.019023 


1.4323 

1.4509 


0.874 

0.019630 

1.4646 

4.5 

0.828 

0.020398 

1.4743 

5.0 

0.788 

0.021313 

1.4805 

5.5 

0.752 

0.022366 

1.4838 

6.0 

* jC 


SECTION  6 


SUMMARY  OF  FORMULAS  FOR  SYSTEMS  ANALYSIS 


Suppose  that  one  wishes  to  calculate  peak  intensity  on  target  for 
a given  set  of  conditions.  The  required  steps  are: 

(1)  Choose  beam  shape  and  set  C£,  a,  m,  m ' , and  m". 


Infinite  Gaussian  0.010590  1.6419 
Truncated  Gaussian  0.028727  1.3715 


Uniform  Beam 


0.014264  1.1777 


0.6366  1 

0.9166  0.8893 
0.9202  1.124 


(2)  Evaluate  nondimensional  numbers. 


Absorption  number  N - 


Slew  number 


Fresnel  number 


•9  P2tNA 

Distortion  number  N_  » (2.333  x 10  y)  — r-— 


(3)  Evaluate  linear  dispersions 


rQ  - 3.017 (k2cjzt)"3/5 

D2  * 8R^,  3.7R2,  or  4R2  for  infinite  Gaussian,  1/e 

truncated  Gaussian,  or  uniform  beam,  respectively 


2 

o _ appropriate  to  hardware 
2 

aw  appropriate  to  hardware 

4 - 4*  4*  4*  4 


(4)  Evaluate  effective  beam  qualities. 


,22,  2 A 2.1/2 

(m  + o,) 
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